Effects of chemical pressure on the Fermi surface and band dispersion in the 

electron-doped high-T c superconductors 
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We have performed angle-resolved photoemission spectroscopy (ARPES) measurements and 
first-principles electronic structure calculations on the electron-doped high-Tc superconductors 
Lni.ssCeo.isCuCU (Ln — Nd, Sm, and Eu). The observed Fermi surface and band dispersion 
show such changes that with decreasing ionic size of Ln 3+ , the curvature of the Fermi surface or 
—t'/t decreases, where t and t are transfer integrals between the nearest-neighbor and next-nearest- 
neighbor Cu sites, respectively. The increase of t with chemical pressure is found to be significant, 
which may explain the apparently inconsistent behavior seen in the hole-doped L^-mSr^CuCU un- 
der epitaxial strain [M. Abrecht et al., Phys. Rev. Lett. 91, 057002 (2003)]. A gap due to the 
antiferromagnetism opens even in the nodal region for the Sm and Eu compounds, and the gap size 
increases in going from Ln — Sm to Eu. 

PACS numbers: 74.72. Jt, 71.20.-b, 79.60.-i, 74.62.Fj 
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Since the discovery of the high-T c superconductors 
(HTSCs) d, a large number of studies have been per- 
formed in order to obtain higher critical temperatures. 
Among them, pressure effects have attracted much at- 
tention because it causes a dramatic increase of T c in 
many systems. For example, in HgBa2Ca2Cu304, which 
has the highest T c among HTSCs, 2T- rises from 135 K 
to 164 K in a hydrostatic pressure [2(]. In addition to 
mechanical pressure 

01 II H H H, the effect of epitax- 
ial strain in thin films grown on single crystalline sub- 
strates [1, Eoj] , and the effect of "chemical pressure" , 
where the lattice constants are varied through substitu- 
tion of ions with different ionic radii [111, EH, EH, have 
been studied so far. The effects of epitaxial strain and 
"chemical pressure" on T c are consistent with the me- 
chanical pressure in many cases [!, 0, IE H, 0, H, S EH , 
but the mechanism of the T c changes has not been un- 
derstood yet. Also, differences in T c between different 
cupratc families are expected to provide crucial informa- 
tion about the mechanism of high-T c su per conductivity 
and have been discussed extensively [iH EH- According 
to first-principles calculations, the shape of the Fermi sur- 
face strongly depends on the distance drju-Oap between 
the copper and the apical oxygen. A long (icu-Oap leads 
to a strong curvature of the Fermi surface, namely, a large 
value of — t'/t [l4[, where t and t' denote transfer integrals 
between the nearest-neighbor and next-nearest-neighbor 
Cu sites, respectively, in the single-band tight-binding 
model. However, the angle-resolved photoemission spec- 



troscopy (ARPES) results of La 2 _ !C Sr a: Cu04 (LSCO) un- 
der epitaxial strain have demonstrated that compressed 
in-plane lattice constant and hence increased (icu-Oap re- 
sulted in a decrease of —t'/ t flrjfTTl ]. contrary to the ma- 
terial dependence of —t'/t [1J]. That is, the relationship 
between the crystal structure and the electronic structure 
of cuprates appears to be more complicated. 

In order to clarify the above issues, we have performed 
ARPES measurements and first-principles electronic 
structure calculations on the electron-doped HTSCs 
Ln2_ x Ce :r Cu04 {Ln = Nd, Sm, and Eu), where with 
decreasing ionic radius of Ln 3+ from Ln = Nd to Eu, 
"chemical pressure" increases, that is, both im^plane and 
out-of-plane lattice constants become small [TJ E3 an d 
T c decreases [ill, Ell ■ Since these materials have no api- 
cal oxygen, the shape of the Fermi surface is not related 
to the <icu-Oap, but determined by other contributions 
such as the in-plane lattice constant. Furthermore, no 
empirical relationship between —t'/t and T c has been 
known for the electron-doped cuprates. It is therefore 
highly desired to investigate the electronic structure of 
the I/n2- x Cea;Cu04 by ARPES. The present results show 
that chemical pressure reduces the curvature of the Fermi 
surface as in the case of the ARPES results on LSCO thin 
films (ltij . Further analysis of the band dispersions has 
revealed that the increase of t as a result of the decrease 
of the in-plane lattice constant has a significant effect on 
the Fermi surface shape. 

High-quality single crystals of optimally doped 
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Ndi.8 5 Ceo.i5Cu0 4 (NCCO), Sm 1 . 85 Ceo.i 5 Cu04 (SCCO) 
and Eui.85Ce .i 5 CuO4 (ECCO) were grown by the trav- 
eling solvent floating zone method. The T c 's of NCCO, 
SCCO, and ECCO were ~22 K, ~16 K, and K, re- 
spectively. There was slight deviation in the Ce content 
and/or the oxygen stoichiometry from the indicated com- 
position, as reflected on small differences in the Fermi 
surface areas among NCCO, SCCO, and ECCO, as we 
shall see below. The ARPES measurements were per- 
formed at bcamlinc 28A of Photon Factory (PF), Insti- 
tute of Materials Structure Science, High Energy Acceler- 
ators Research Organization (KEK), using incident pho- 
tons with energy of 55 eV. We used a SCIENTA SES-2002 
electron-energy analyzer in the angle mode, and used a 
five-axis manipulator [l8j |. The total energy resolution 
and momentum resolution were 15 meV and 0.017T, re- 
spectively. Samples were cleaved in situ under an ul- 
trahigh vacuum of 10~ n Torr to obtain clean surfaces, 
and were measured at ~ 10 K. The Fermi edge of gold 
was used to determine the Fermi level (Ep) position and 
the instrumental resolution before and after the ARPES 
measurements . 

Figure 1 shows plots of ARPES intensity at Ep in 
NCCO, SCCO, and ECCO in two-dimensional momen- 
tum space. The suppression of the intensity is seen near 
the "hot spots", i.e., the intersecting points of the para- 
magnetic Fermi surface and the antiferromagnetic Bril- 
louin zone boundary [see the inset in Fig. 1(c)]. This 
suppression is due to the scattering of electrons at E-p 
by antiferromagnetic fluctuations or by (quasi)static an- 
tiferromagnetic correlation with wave vector (jr, it) as dis- 
cussed in previous ARPES studies [HH^. We present 
energy distribution curves (EDCs) in the nodal direction 
in Fig. 2, and ARPES intensity plot in E — k space in the 
nodal direction in Fig. 3(a), (b), and (c). In Fig. 2, the 
gap or leading edge (LE) shift is observed for SCCO (LE 
shift: - 5 meV) and ECCO (LE shift: - 30 meV) but 
not for NCCO, consistent with a recent ARPES study on 
Sm 1 . 86 Ce . 14 CuO4 M- 

White circles in Fig. 1 represent the peak positions 
of the momentum distribution curves (MDCs) at E-p and 
represent the Fermi surface or remnant Fermi surface. As 
seen from Fig. 1, the curvature of the Fermi surface in 
NCCO is the strongest among the three materials. We 
quantitatively evaluated the difference of the curvature 
among the three compounds using a tight-binding model 
as follows. We used two-dimensional antiferromagnetic 
tight-binding model: 



e — fi = Eo ± y AE 2 + 4£ 2 (cos k x a + cos k y a) 2 
— 4t'cos fc x acos k y a — 2t" (cos2fc x a + cos2k y a), 

where t, t', and t" are transfer integrals between the 
nearest-neighbor, second- nearest-neighbor, and third- 
nearest-neighbor Cu sites, respectively, eo represents the 
center of the band relative to the chemical potential /z, 
and 2AE is the potential energy difference between the 
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FIG. 1: (Color online) ARPES intensity within ±30 meV 
of the Fermi level (-Ef) plotted in momentum space for 
Nd 1 . 85 Ceo.i 5 Cu0 4 (NCCO) (a), Sm 1 . 85 Ceo.i5Cu0 4 (SCCO) 
(b), and Eui.85Ceo.i5Cu04 (ECCO) (c). The data were taken 
over a Brillouin zone octant and symmetrized with respect to 
the (0,0)-(n,n) line. White circles show the peak positions of 
momentum distribution curves (MDCs) at Ef, indicating the 
shape of the Fermi surface or underlying Fermi surface. Solid 
red curves and dashed pink curves show the Fermi surface 
obtained by tight-binding fit to the ARPES data assuming 
the paramagnetic and antiferromagnetic band structures, re- 
spectively. Inset: Schematic diagram of the hot spot. Black 
curve and red dashed line represent the Fermi surface and the 
antiferromagnetic Brillouin zone boundary, respectively. 



spin-up and spin-down sublatticcs. 

First, we fitted the calculated Fermi surface to the set 
of white circles in Fig. 1 by adjusting the parameters 
—t'/t and —£o/t. We note that the curvature of the Fermi 
surface depends on the ratio —t'/t and not on each of t 
and t' , and only very weakly on —So/t. Assuming that 
-t" It' = 0.50, we obtained -t'/t = 0.20, 0.12, and 0.11 
for NCCO, SCCO, and ECCO, respectively, and e /t = 
-0.12, 0, and -0.05 for NCCO, SCCO, and ECCO, re- 
spectively. The Fermi surface area was allowed to de- 
viate slightly from what is expected for 15 % Ce dop- 
ing because of non-stoichiometry. Note that when we 
performed the same analysis retaining only the second- 
nearest neighbor Cu hopping (t" = 0), we obtained 
-t'/t = 0.40, 0.23, and 0.21 for NCCO, SCCO, and 
ECCO, respectively, showing a similar tendency to the 
case of —t"/t' — 0.50. The present fitted results demon- 
strate that the curvature of the Fermi surface monotoni- 
cally decreases from NCCO to SCCO to ECCO. 

Next, we evaluated the absolute values oft and t' by fit- 
ting the calculated band dispersion to the MDC and/or 
EDC peak positions with —t'/t and —So/t fixed at the 
above values. Because the position of the "flat band" at 
~ (tt, 0) is sensitive to —t', it can be used to determine 
t' and t separately. Figure 3 shows the E — k plot of 
ARPES intensity for NCCO, SCCO, and ECCO along 
two different cuts as shown in the inset. Panels (d)-(f) 
show that the positions of the band around the (tt, 0) 
point becomes shallower in going from NCCO to ECCO, 
indicating that — t' decreases. For the plot of the antifer- 
romagnetic tight-binding bands, we have assumed that 
the spectral intensity denoted by the size of the point 
is proportional to the projected weight of the antiferro- 
magnetic band to the paramagnetic band. In the nodal 
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FIG. 2: (Color online) Energy distribution curves (EDCs) for 
NCCO (a), SCCO (b), and ECCO (c) around the nodal k F 
point. Blue thick EDCs represent EDCs at fcp. The direction 
of the cut is shown in the insets, (d) EDCs of NCCO, SCCO 
and ECCO at the nodal k F point. 



region [panels (a)-(c)] the experimental band dispersion 
as well as spectral weight is well explained by the antifer- 
romagnetic band while around the (ir, 0) region [panels 
(d)-(f)], the experimental one cannot be explained by the 
antiferromagnetic band well but rather by the paramag- 
netic band (AE = 0). These fitting results may indicate 
that the fc-dependence of the antiferromagnetic gap AE 
exists beyond standard band theory of the antiferromag- 
netic state, as predicted by a recent variational Monte- 
Carlo study H3- 

In order to see the effect of "chemical pressure" on the 
band structure from first-principles calculation, we have 
carried out muffin-tin-orbital (MTO) based NMTO cal- 
culations [23| within the framework of local density ap- 
proximation (LDA) and extracted the transfer integrals 
corresponding to t, t' , and t" in the tight-binding model. 
For this purpose, an effective Cu d x 2_ y 2 basis was defined 
by means of the downfolding procedure, by integrating 
out all the degrees of freedom related to Ln, O and Cu 
except for Cu d x 2_ y 2. The effective basis, constructed 
in this manner, serves the purpose of the Wannier-like 
function corresponding to the single band crossing the 
E-p. The real-space Hamiltonian defined in the basis of 
these effective, Wannier-like orbitals provides the infor- 
mation about the various transfer integrals, connecting 
various Cu sites. This method has been applied success- 
fully in the case of hole-doped cuprate compounds [l4|. 
Calculations have been carried out for Ln =Nd and Sm. 
The calculated -t'/t turned out to be 0.34 (t = 0.41 eV, 
t' = -0.14 eV) for NCCO and 0.29 (t = 0.44 eV, t' = -0.13 
eV) for SCCO. Such moderate change in —t'/t, compared 
with the hole doped cuprate family is expected consider- 
ing the absence of apical oxygen in the T' structure, one 
of the key controlling factor in determining —t'/t [l4| . 
Though the quantitative estimates of —t'/t are somewhat 
different from those obtained from experimental data, the 
trend is correctly described as shown in Fi g. 4. Accord- 
ing to cellular dynamical mean-field theory [2J|, —t'/t is 
reduced by electron correlation over the LDA estimate 
for electron-doped cuprates and increased for hole-doped 
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FIG. 3: (Color online) ARPES intensity plot in energy- 
momentum space in the nodal and antinodal regions of the 
Brillouin zone compared with tight-binding energy bands, 
(a), (d): NCCO, (b), (e): SCCO, and (c), (f): ECCO. The 
cuts are shown in the inset. The spectral intensity of the 
antiferromagnetic band is represented by the size of squares. 



systems, with the effect being strong for electron-doped 
case at low doping level. This may explain the dispar- 
ity in the quantitative values of —t'/t as seen in Fig. 4, 
where for LSCO a good agreement in —t'/t between the 
LDA and the experimental estimate was observed. 

Considering the nearest neighbor Cu-Cu transfer inte- 
gral t to be given by the relationship t — 2(t pc [) 2 / {£d— £ P ), 
where t p d and the Cu 3d - O 2p transfer in- 

tegral and the onsite energy differences, respectively, de- 
fined within the three band model, the increase of t in 
going from Nd to Sm is contributed both by the increase 
of tpd due to contraction of the lattice, as well as due 
to decrease of (e<j — e p ) caused by the replacement of 
Nd by Sm. Our NMTO-downfolding calculations keep- 
ing oxygen-per degrees of freedom active in addition to 
Cu d x 2_ y 2 gives a 8% change in t p d, while the rest is con- 
tributed by the change in (e^ — e p ). The decrease of the 
second nearest neighbor hopping, t' in moving from Nd 
to Sm is, on the other hand, contributed by the decrease 
in the 0-0 hopping, t pp as obtained in our three band 
calculations, presumably caused by the different cation 
covalency effect between Nd and Sm. 

We now turn our attention to the case of strained 
LSCO thin films. In LSCO, the reduced in-plane lat- 
tice constant by compressive strain lead to the decrease 
of —t'/t [16[, as in the present electron-doped system. 
Hence, the experimental results indicate that the reduced 
in-plane lattice constant resulted in the reduced value of 
—t'/t both under chemical pressure and epitaxial strain 
and both for the hole-doped and electron-doped com- 
pounds. Also, since the electron-doped HTSCs have no 
apical oxygen, the shape of the Fermi surface in our ex- 
periment is not related to (icu-Oap but is related to the 
in-plane lattice constant including the effect of 0-0 hop- 
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FIG. 4: (Color online) (a), (b): Tight-binding parameters 
as functions of in-plane lattice constant. For definition, see 
text. As the in-plane lattice constant decreases, — tl and —t'/t 
decreases while t and AE increase. Dashed lines are results 
from the LDA calculations. LSCO data are taken from Refs. 

MM. 



ping. Therefore, we consider that in the strained LSCO 
film the effects of in-plane lattice constant is stronger 
than those of dcu-Oap- 

Finally, we discuss the relationship between —t'/t val- 
ues and the superconductivity. In the hole-doped HTSCs, 
there is an empirical trend that the larger —t'/t is, the 
higher T c becomes [14]. The present results on —t'/t 
in the electron-doped system follow the same trend, but 
considering LDA estimates of —t'/t, the predicted T c 
would be higher, ~ 90 K and - 70 K for NCCO and 
SCCO, respectively. If we employ the experimental val- 
ues for —t'/t, instead, then T c would be ~ 40 K, 25 K for 
NCCO and SCCO, in somewhat better agreement with 
the empirical trend, but still overestimating the differ- 
ence in T c between NCCO and SCCO. The gap open- 
ing in the nodal direction may explain the generally low 



T c 's in the electron-doped systems compared with hole- 
doped ones as well as the systematic suppression of T c 
with chemical pressure because electrons near the nodal 
point significantly contribute to the superconductivity, 
according to the Raman scattering studies [26| and varia- 
tional cluster calculations [27[- In addition, in going from 
NCCO to ECCO, the compression of the out-of-plane 
lattice constant becomes strong compared with that of 
the a-axis 11] . The decrease of the inter-layer distance 
may cause the development of three-dimensional antifer- 
romagnetic order, leading to the large AE. Further stud- 
ies are necessary to clarify whether the effect of variation 
in —t'/t or that in AE is more important for the change 
of T c . 

In conclusion, we have performed ARPES and first- 
principles electronic structure calculation studies of 
NCCO, SCCO, and ECCO in order to elucidate the vari- 
ation of electronic structure by "chemical pressure" . As 
the in-plane lattice constant decreases, —t'/t decreases 
and t increases, consistent with the previous ARPES re- 
sults on strained LSCO films. This suggests that the 
variation of the in-plane lattice constant has a great in- 
fluence on the electronic structures. We consider that 
the decrease of T c with increasing "chemical pressure" is 
attributed to the change in the Fermi surface and/or the 
gap opening around the nodal point. 
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